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Abstract 

We study the Landau level localization and scaling properties of a disordered two- 
dimensional electron gas in the presence of a strong external magnetic field. The 
impurities are treated as random distributed scattering centers with parameterized po- 
tentials. Using a transfer matrix for a finite-width strip geometry, we calculate the 
localization length as a function of system size and electron energy. The finite-size 
localization length is determined by calculating the Lyapunov exponents of the trans- 
fer matrix. A detailed finite-size scaling analysis is used to study the critical behavior 
near the center of the Landau bands. The infiuence of varying the impurity concentra- 
tion, the scattering potential range and its nature, and the Landau level index on the 
scaling behavior and on the critical exponent is systematically investigated. Particular 
emphasis is put on studying the effects of finite range of the disorder potential and 
Landau level coupling on the quantum localization behavior. Our numerical results, 
which are carried out on systems much larger than those studied before, indicate that 
pure (5-function disorder in the absence of any Landau level coupling gives rise to non- 
universal localization properties with the critical exponents in the lowest two Landau 



levels being substantially different. Inclusion of a finite potential range and/or Landau 
level mixing may be essential in producing universality in the localization. 

PACS numbers: 73.40.Hm,71.30.+h,73.50.Jt 
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I. Introduction 



The quantum Hall effect [1] phenomenon is now well accepted to be closely related to the 
strong-field Landau level localization problem in a disordered two-dimensional(2D) electron 
system. In a disordered 2D system, without any external magnetic fields, it is well-known 
that all electronic states are localized [2] and strictly there is no metal-insulator transition. 
In the presence of a strong external magnetic field the energy spectrum of a 2D system is a set 
of impurity-broadened Landau bands. While the electronic states in the tails of the Landau 
bands remain localized, the states at Landau band centers delocalize and become extended 
with a mobihty edge Ec near the center of each Landau band. When the Fermi energy is 
in the localized regime, the Hall resistance (p^^y) is quantized with vanishing dissipative or 
longitudinal Yesistance{pxx = 0) at T = 0. When the Fermi energy moves to the delocalized 
regime, there is a metal-insulator transition at Ec, and the system behaves as a metal with 
non-zero p^x and non-quantized pxy which takes us from one quantized Hall plateau to the 
next. This transition regime, where p^x is non-zero and Pxy is jumping from one plateau 
to the next, plays an important fundamental role in localization and quantum Hall studies. 
Experimentally [3-9], this transition regime becomes sharper and narrower as T decreases 
with the T = limit thought to be infinitely sharp with only the state at E = Ec being 
extended and all other states localized. Existence of these extended states(a set of measure 
zero at T = 0) near the Landau level centers is of crucial importance to the quantum 
Hall effect phenomenon and, conversely, the quantum Hall effect demonstrates the existence 
of extended states (and, of metal-insulator transitions) in the strong-field 2D system. In 
this paper, we provide a detailed numerical study of the localization properties of the two 
lowest Landau levels in the strong-field 2D system using a finite-size scaling analysis. There 
exist several recent studies of this type in the hterature [10-13]. Utihzing the advantages 
provided by massively parallel computing machines, we have been able to go to system sizes 
substantially larger (by a factor of four) than those hitherto existing in the literature. In 
addition, our emphasis, in contrast to the earlier numerical studies [10,11] which mostly 
concentrate on the lowest Landau level in the presence of zero-range disorder potential, is 
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on the effects of finite range disorder potentials and Landau level coupling effects on the 
localization properties of the two lowest Landau levels. We are particularly interested in the 
issue of universality {i.e. whether the critical localization exponents are always the same 
independent of Landau level index, type of disorder potential, etc.) in strong-field Landau 
level localization , a topic which has created some controversy both in the experimental and 
in the theoretical literature on the subject [3-14]. 

The critical localization behavior of the quantum Hall effect can be summarized as follows. 
The quantum Hall resistance p^y and the longitudinal dissipative resistance p^x of a strong- 
field 2D disordered system are measured at different temperatures(T) in the transition regime 
by varying the magnetic field(5). In the transition regime it is found [3,4] that the maxima of 
dP-pxy/dB'^ diverge as T~'^'^{n — 1, 2, 3) and the width of p^x peak, AS, is proportional to T'^. 
Here, AB is the width of the transition regime where pxx is non-zero and pxy is unquantized. 
The exponent k of the temperature dependence arises from a competition between two 
microscopic length scales at a finite T [14]. One is the inelastic scattering or the phase 
breaking length Lj„ which scales with temperature with an exponent p(Lj„(T) oc T~^/^). 
The other is the localization length ^ for the electronic state at the Fermi energy, which 
scales with energy with an exponent v. The localization length diverges as the Fermi level 
approaches the mobility edge critical energy Ef. in the middle of each Landau band i.e. 
^{E) oc| E — Ec \~'^. When Ep is at the tails of Landau bands, ^ <S Lj„, the system 
is insulating and localization plays the major role in producing quantum Hall plateaus. 
When Ep is near the critical energy, ^ ^ Li„, inelastic scattering length acts as a cut-off 
making the system behave metallic. Using scaling arguments [14] it is easy to combine these 
two exponents to show that the temperature dependence of pxx and pxy in the transition 
regime(ie. in the metallic phase with Ep ~ Ec) would be controlled by a composite exponent 
K — p/2i/, which is directly measured experimentally by studying the Pxx peak and the 
associated dpxy/dB. (We have, in fact, recently shown [15] that the scaling analysis can 
be extended even to the pxx minima, ie. where pxy is quantized, provided the conduction 
mechanism is activated transport and not variable range hopping.) A key question in this 
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problem is whether the exponent n is universal or not. Earlier experimental investigation 
[3,4] involving studies of the temperature dependence of p^x peak and the associated dp^y/dB 
was controversial with k being independent of Landau level index in one study and dependent 
(by a factor of two) on the Landau level index in the other. There are also some reports [5] 
of the dependence of the localization exponent on the impurity scattering strength. In our 
opinion, the most thorough and reliable temperature dependent experiments tend to indicate 
a universal (Landau level independent) k. The problem with experimental measurements of 
K is that it is a composite scaling exponent, k — p/2i/, being dependent on the localization 
exponent p and the phase breaking exponent p. While calculations, such as the one carried 
out in this paper, can provide a reasonably accurate numerical value of i/(for non-interacting 
electrons), the value of the inelastic scattering exponent p in the strong- field situations is 
simply not known. Comparison between theory and temperature-dependent localization 
experiments has been carried out mainly under the assumption that p ^ 2 which is the 
corresponding zero-field clean limit result for electron-electron scattering in two dimensional 
Fermi hquids [10]. (We mention that under the assumption p — 2, one gets k, 0.4 using 
the numerically calculated vo ~ 2.3, obtaining excellent agreement between theory and 
experiment.) There are, however, some difficulties with this "standard" theoretical analysis 
of the localization data. First, we do not see any justification for applying the zero-filed 
theoretical result for p to the strong-field situation. Second, there is no justification for using 
the clean limit result (p ~ 2) rather than the disordered 2D result (p ^ 1) except perhaps 
the(fortuitous?) "agreement" obtained between theory and experiment using p — 2. Third, 
even if v is universal(ie. Landau level independent, etc.), it is very hard to see how p could 
be universal - for example, in the zero-field situation p = 1(2) depending on whether the 2D 
system is dirty (clean). Because universality of k requires universality of both p and u, we 
believe that a complete understanding of the experimental results awaits a theory for inelastic 
scattering in the strong-field 2D system. In this respect, we note that temperature dependent 
experiments measuring the width AB of Pxx peaks and of the associated dp^y/ dB show that k 
in GaAs samples[9] achieve the "universal" value (~ 0.4) only for T < 200mK whereas earlier 
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similar experiments [3] on InP heterostructures produce the universal k ~ 0.4 for T ~ AK. 
Thus, the size of the critical scaling regime (which is determined by an interplay between 
localization and inelastic scattering) is clearly sample-dependent. There are also systematic 
studies of Si-MOSFETs[4,6,7] and GaAs heterostructures [5] in the literature which find that 
the measured value of k is non-universal, with the reported dependence of k on both the 
Landau level index and sample quality. A recent size-dependent experimental localization 
study[8] of narrow GaAs heterostructures, which claims separate experimental measurements 
of 1/ and p, concludes that while 2.3) is universal (independent of both Landau level 
index and the sample quality), p varies from 2.7 to 3.4, depending on the Landau level index 
and the sample quality. Thus, k itself («; = p/2u) may not be universal. We emphasize 
that our work is solely on an accurate numerical calculation of the localization exponent v 
and the universality being discussed in this paper refers entirely to the localization aspect 
of the problem(z.e. we are concerned with whether or not v is independent of Landau 
level index, the nature of the impurity scattering potential, etc.). We have nothing to say 
about the exponent p which is often needed in understanding the experimental quantum 
Hall localization data. 

The interesting recent experimental attempt [8] , alluded to above, circumvent the problem 
of having a composite exponent k by direct and independent measurement of i/(and p) using 
experimental 2D strip systems of various widths warrants a discussion. The competing length 
scales in this type of experimental "finite-size scaling" studies are the localization length 
and the system size (with the inelastic scattering length, which is not playing an active role, 
assumed to be a constant, larger than the system size, because the temperature is held fixed 
in these experiments) . Such direct measurements of v show that it is indeed universal with a 
value u a; 2.3 in excellent agreement with the existing (and our own) theoretical calculations. 
While the basic premise of these experiments is novel and this procedure certainly eliminates 
the difficulty associated with k{— p/2v) being a composite exponent, we feel that there are 
some inherent problems in the interpretation of these experimental results as well. Suppose 
the temperature is such that Lg > where Lg is the finite size of the experimental 
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system. In that situation, obviously Lg is irrelevant for the localization behavior of the 
system and the competition between the localization length ^ and the inelastic scattering 
length Lin controls the scaling behavior. Thus, for the finite-size scaling experiments to 
be meaningful, one must have Lj„ >> Lg, which is precisely the mesoscopic conductance 
fluctuation regime [16] where the system is manifestly non-self-averaging and one should see 
quantum interference induced conductance fluctuations because the system size is smaller 
than the phase coherence length. Such strong mesoscopic fluctuations are not present in 
the experimental results and, therefore, it is unclear to what extent the necessary condition 
Lin » Lg is satisfied. We emphasize that, while the weak-field conductance fluctuation 
phenomenon is currently well understood, there has been very little work on strong-field 
conductance fluctuation phenomenon, and, in principle, it is possible that there is a very good 
(but presently unknown) reason why conductance fluctuations are absent in the strong-fleld 
2D systems even in the mesoscopic L^n » Lg regime. Without such a theory explaining 
the absence of strong-fleld mesoscopic conductance fluctuations, we feel that a complete 
understanding of the novel flnite-size scaling experiments remains elusive. 

In this paper, we study localization and scaling properties of a disordered two-dimensional 
system of flnite width in the presence of a strong magnetic fleld. The impurities are treated 
as random distributed scattering potential centers. By choosing an appropriate cut-off in 
the scattering range, we can divide our system into cells with only the nearest-neighbor 
intercell scattering being appreciable. This allows us to set up an interccll transfer matrix 
for evaluating the exact single particle non-interacting electron wavefunction in the presence 
of disorder. The flnite-size localization length of the electron is determined by calculating 
the Lyapunov exponents of the transfer matrix on a massive parallel processing machine. 
Finite-size scaling method is applied to study the critical behavior near the center of the 
Landau bands. The influence of varying the impurity concentration, the scattering potential 
and the Landau level index on the critical exponent is systematically investigated. A short 
report of our work has earlier appeared in the literature[17]. 

The rest of the paper is organized in the following way: in Sec. II we give the details of 
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our theoretical method with the relevant formulas and equations; in Sec. Ill we present our 
numerical results; we provide a discussion with the conclusion in Sec. IV. 

II. Theoretical Method 

For a non-interacting 2D electron gas in the presence of a strong external magnetic field 
B = (0, 0, B) perpendicular to the 2D x — y plane, the single electron energy spectrum is 
given by a set of discrete Landau levels (denoted by the index N — 0,1,2, ...): 

Er,^iN+^)hu;,, (1) 

where cUc — eB/mc is the cyclotron frequency, and m the electron effective mass. Each 
Landau level is highly degenerate with a macroscopic degeneracy given by (27rZ^)~^ per unit 
area where Ic = {Uc/eBy/'^ is the magnetic length. By choosing the Landau gauge for the 
vector potential A = (0, 5a;, 0), and periodic boundary conditions in the y direction, the 
electron wavef unction becomes a simple harmonic oscillator eigenf unction in the x direction, 
and a free-electron plane wave in the y direction. The one electron Landau wavefunction of 
the pure system can then be written as 

^»-^-*"-"'^)-*'^'- 

Here M is the width of the 2D system in y direction, if„ is the Hermite polynomial, and 
k — "^n where n is an integer. Note that the simple harmonic oscillator has a displaced 
center with the oscillation center sX x = X = kll- We choose these one electron Landau 
eigenstates as a complete basis in the Hilbert space for our localization study of the disordered 
2D system. 

We consider a two-dimensional disordered system of length L (along the x direction) and 
width M in the strip geometry with periodic boundary conditions applied in the y direction. 
We use lo — \/2t^Ic as the length unit in this problem. The single electron Hamiltonian in 
the presence of disorder can be written as 

H = Y,\NX > {N + )-)huc <NX\ \ NX ><NX\V \ N'X' >< N'X' \ (3) 

NX ^ NX N'X' 
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where | NX > is the Landau state for the A^th Landau level with the center of oscillation 
X taking discrete values with spacing 1/M due to the periodic boundary conditions. The 
disorder potential V, which is treated as the one electron potential arising from a random 
distribution of impurity scattering centers, is given by 

V{r)^Y.V^i^-^^) (4) 

i 

where is the location of the ith impurity scattering center and Vi depends on the type of 
scattering potential. We consider two models for the random disorder potential, namely, 
the short-range potential , Vi — VoS{r — rj), and, the finite range Gaussian potential, 
Vi = ;^e-l'-""'l . The potentials are chosen to be all attractive, or all repulsive, or an 
equal mixture of both. The concentration of impurities, which along with Vq determines the 
scattering strength, is denoted by q. An important numerical quantity in our calculation is 
Icutoff which determines the length scale above which the impurity scattering potential does 
not connect the different Landau states, i.e. we only need to consider scattering between 
states X,X' satisfying | X — X' \ < Icutoff- Obviously Icutoff is determined by the characteris- 
tics of the Landau wavefunctions and the impurity potential and concentration. Our criterion 
to determine Icutoff is that the matrix element < NX \ V \ N'X' > ior \ X — X' \> Icutoff 
should be much smaller (less than 1%) than that iox \ X — X' \< Icutoff- In our calcula- 
tions, we choose Icutoff — 1 for the N — Landau level, and Icutoff — 2 for higher Landau 
levels (A^ — 1). Our choice of I cutoff is consistent with the others existing in the literature 
[10,11]. In addition to the inter-Landau-level energy separation hujc, there is another energy 
scale in this problem which is the Landau level broadening Fjv determined by the strength 
of the impurity-scattering potential. Within the self-consistent-Born-approximation(SCBA), 
and in the strong field limit {i.e. no Landau level coupling), the density of states(DOS) for 
the short-range impurity potential is given by [18]: 



D{E) = 2j:{nT)-' 



1 - 



E-E 



N 



1/2 



(5) 



N=0 

where F = 2Vo^/cl <^ hUc- Defining 7 = V /hwc, we note that 7 << 1 defines the weak 
disorder strong-field limit when scattering to other Landau levels may be ignored. In some 
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situations we would consider the strong disorder case, 7 ~ 1, where we have to take into 
account Landau level coupling in our calculations. We choose hcUc as the unit of energy in 
our calculations (which is the most obvious natural choice). 

Our localization calculation proceeds in two steps: First, we obtain the localization length 
for a 2D strip of a finite width M and an essentially infinite length of L = 10^(>> M — 
4 — 256) by calculating the system Lyapunov exponent on a massively parallel processing 
machine; then, we carry out a finite size scaling analysis to obtain the localization length 
in the M — > 00 limit. In order to calculate the localization length for a 2D strip of length 
L and width M, we start by dividing the system into cells with length Icutoff and width M 
such that the next-nearest-neighbor intercell interaction can be neglected. The number of 
states in each cell for each Landau level is n = Icutoff M. Expanding the wavefunction |^ > 
of the disordered system for energy E in the complete basis of Landau wavef unctions, we 
get 



with 



* >= ^ aNi I NXi > 

m 



H \ >^ E \ > . 



(6) 



(7) 



Combining Eqs.(l),(3),(6) and (7) we get 



1, 



2n 



{N + -)huJcaN,i+n + N^i+n \ V \ N' X^+k > aN',i+k = EaN,i+n- 

^ N' k=0 



(8) 



Defining column vectors 



and 



an 



(9) 



( A ^ 



ii) 



(10) 
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where Nl is the number of Landau levels being included in our numerical calculation, Eq. (||) 
can be transformed to the following set of linear equations: 

2n 



where 



i+k 



\TluJc 

\hujc 



k=0 






^kn 



(11) 



V ... {Nl - '^)hu, j 
( < 0X,+„ I V I ox,_,, > 



(12) 



+ 



< OXi+„ \ V\Nt^- 1,X;+, > 



\ < NL-l,Xi+n\V \OXi+k> ... <NL-l,X,+n\V\NL-l,X,+k> J 

Equation (pi]) shows that Ai+2n is determined by Ai, Ai+i...Ai+2n-i i-e. we can set up the 
transfer matrix T^^^ such that 

^{i+i) = (13) 
Here the transfer matrix T^^^ is a (2ri x Ni) x (2ri x X^,) dimensional matrix: 



rp{i) 



/ ... 



V 















(14) 



where / is Nl x Nl unit matrix. Note that if Landau level coupling is neglected in the 
calculations (as has been done in all the existing calculations), Nl = 1. We have done 
calculations for single Landau levels [N = 0, 1) with Nl = 1 and for coupling between the 
lowest two Landau levels with Nl = 2. T^*^ is calculated by summing up the contribution 
from each impurity scatterer, and, therefore, depends on the details of the local distribution of 
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the random impurities. For example, for the short-range impurity potential Vi — VoS{r — Vi), 
we obtain the following non-zero matrix elements: 



< OXn, I Vi I IXn, > 



M 



< IXn, I Vi I OXn, > 



< OX^^ I Vi I 0X„, > = 



M 

xe~''[(^'~"i/-'^)^+(^'~"2/M)2]^ ^^^^ 

^^M^i - ni/M){Xi - n2/M)e'2^(ni-n2)2/i/M 
y^Q-Tvl{xi-m/M)^+{xi-n2/M)^] ^ ^-j^g-j 

We do not show here the explicit matrix elements for the Gaussian impurity potential which 
can also be calculated analytically. 

As mentioned before, once the transfer matrix T^*) is formed, we calculate the localization 
length by evaluating the Lyapunov exponent for the transfer matrix. Calculation of the 
Lyapunov exponent for maps such as the one defined in Eq.(13) is standard in non-linear 
dynamics, and has been successfully used in localization calculations in other contexts [19]. 
For the transformation T which transforms a vector ip according to Eq.(13), the Lyapunov 
exponents are given by 

1 I ^(^+1) I 1 I nil tw^^^^ I 

Ly = lim - In-— = lim - In ' , , . (19) 

For our (L, M) strip geometry, with L very large and M small, | tp^^^ |~| ■0(0) |, | tp'^'''^^^ |~| 
ip{L) I, A; = = LM. For localized states, the localization length Am is defined to be 
the characteristic length controlling exponential decay of the wavefunction, i.e. \ ip{L) \ / \ 
ip{0) |~ e-^/^'^, so that 

^" = -Ml,- (^°) 

Mathematically there exist 2n x Nl Lyapunov exponents for the transfer matrix T, but the 
relevant one for the localization length is the smallest negative Lyapunov exponent. 
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We employ the following technique to calculate the smallest negative Lyapunov exponent. 
Denote S as an upper triangular matrix, i.e. Si^j = 0, Q as a lower triangular matrix, i.e. 
Qi<j = with detg = 1. T(i) can be expressed as T^^^ = Q^^^S^^\ then r(2)g(i) = ^(2)5(2) 
and so on. Finally, 

k k 

f[ T« = Q« f[ (21) 



i-l i=l 

The relevant Lyapunov exponent we need is 

k 

In I 

fe— ►oo k 



^ J™ I 11'^^'^ |nxiVz,+l,nxiVz,+l (22) 

1=1 

The problem of computing the Lyapunov exponent now reduces to a simple matrix mul- 
tiplication problem. The problem is, therefore, ideally suited for the parallel processing 
capability of massive parallel processing computers. We have used both CM2 and CMS 
parallel processing computers of the Thinking Machine Co. to carry out our calculations. 
Note that, since scattering from various impurities is independent of each other, we can use 
a parallel processing algorithm to calculate their contribution to T as well. 

Once the localization length, Xm{E), at a particular energy E and for a specific finite 
value of M is obtained from the Lyapunov exponent calculation, we carry out a finite-size 
scahng analysis to go to the 2D hmit(ie. the M — > oo hmit). Finite-size scahng analysis has 
been extensively used in localization calculations and for details we refer to the literature 
[10,11]. The idea is that one should see one-paramctcr scaling for large enough sample width 
M > Mgc and for energies in the critical regime | E — Ec | < £^sc(iiote that the system length 
L is very large). 

Here f{x) is the universal scahng function for the locahzation problem, with ^{E) the lo- 
calization length for the infinite system. This means that in the critical regime a change 
in energy is equivalent to a change in the system size. As we go to large system sizes, the 
localized states far away from the critical region have ^{E) <^ M, and Am ~ so the scahng 
function has the asymptotic behavior of f{x) ~ l/x{x » 1). On the other hand, for the 
states very close to the critical energy E ~ E^ Xm ~ M, and ^ 3> M, so that f{x) ~constant 
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for X <^ 1. By plotting the numerically calculated Xm{E)/M against M on a log-log plot, 
we can calculate ^(E) by using a least-square optimization fit to make all the data points 
collapse on a single "smooth" scaling curve for Xm{E)/M against M/^{E) obeying Eq.(23) 
and the correct asymptotic properties. This finite-size scaling analysis immediately gives 
us i{E), the locahzation length for the infinite system. Using log-log plot of the calculated 
^{E) against E, we can obtain the critical exponent i.e. 

C{E)(x\E-E,\-'' . (24) 

Obeying of both Eq.(23) and (24) together puts rather stringent requirements on the finite- 
size scaling analysis ensuring the correct evaluation of the critical exponent v from the 
optimization procedure. Note that the long strip-geometry topology of our model ensures 
self- averaging for the Lyapunov exponent calculation of the transfer matrix and no additional 
averaging is required. 

III. Numerical Results 

In this section we present detailed numerical results obtained by applying the theoretical 
method described in II. We have carried out calculations on system sizes up to M = 256, 
which is larger than those studied before [10,11]. The sample length L(~ 10^) is chosen 
such that the numerical error in the calculation of Lyapunov exponent is less than 1%. We 
consider a wide range of impurity concentrations from = 2 to = 32. We consider 
the short-range (5-function scattering potential and Gaussian long-range potentials with the 
potential range d/lc — 0.5 and 1. For the weak disorder case, we carry out calculations for 
single(ie. Nl — 1) Landau levels N — 0,1, and for the strong disorder case, we take into 
account Landau level coupling(iV2, = 2) between N = and 1 Landau levels. 

In Figs.l and 2 we present our numerical results for Landau levels A'" = and 1 respec- 
tively in the weak disorder (ie. neglecting Landau level coupling), symmetric(ie. with equal 
number of repulsive and attractive 5-function scatterers) case. The weak disorder results 
for different impurity concentrations with symmetric 5-function scatterers are summarized 
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in Table 1. The calculated localization critical exponent v depends strongly on the Landau 
level index, but varies little with the impurity concentration. This finding of a strong varia- 
tion of V with the Landau level index N is consistent with earlier findings in the literature 
[10,ll](using substantially smaller system sizes). While our numerical results show reason- 
ably good one-parameter scahng (Fig. 1(b) and Fig.2(b)) properties for both the Landau 
levels, the asymptotic value of the scaling function seems to depend on the Landau level 
index. In our scaling plots, the saturation values for M/^ <^ 1 are quite different for the 
different Landau levels (for = 0, Xm/M 1; for = 1, Xm/M ~ 2). The insensitivity of 
the scaling function to a variation in the impurity concentration Cj is, however, quite robust. 

For random distributed weak-disorder repulsive 6 impurities, the calculated critical expo- 
nents are presented in Table 2 and a typical set of scaling data is shown in Fig.3. In this 
asymmetric situation with all repulsive scatterers, the critical energy Ec is shifted away from 
the exact center of the Landau band, but as the impurities get denser, the critical energy 
approaches the Landau band center. In the dilute impurity limit, the critical exponents are 
asymmetric, being somewhat different for the lower and the higher energy branches of the 
Landau band. It is known that in the low impurity concentration case, the density of states 
(DOS) is asymmetric with the maxima of DOS shifting away from the center of the Lan- 
dau band (for repulsive scatterers, the maxima shift to lower energies) [20]. We believe the 
asymmetry in our calculated critical exponents arises from this finite impurity concentration 
effect, and should disappear if Landau level coupling is included in the calculation. 

Our results for the values of u in the presence of Gaussian impurity scattering potential 
(without any Landau level coupling) is summarized in Table 3, and some typical scaling 
plots are shown in Fig.4 (for N = 0) and Fig.5 (for N = 1). Basic trends are similar to those 
found in the 5-function case with one important difference. While for the N = Landau 
level, the critical exponent remains essentially a constant {u ~ 2.2) as the potential range is 
varied, the exponent for the N — 1 Landau level clearly decreases with increasing the range 
of the impurity potential. In particular, the value of u decreases from around 5.5 for zero- 
range potential to around 2.5 for finite range potentials (the value of i/q, however, remains 
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constant around 2.2). Based on this rather hmited numerical data, one may speculate that 
universality is restored by finite-range scatterers and that 5-function scattering is quite a 
special case. 

In the presence of strong disorder, the Landau level broadening F is comparable to inter- 
Landau-level energy difference hcUc, and, therefore, we must include Landau level coupling in 
the calculation. We consider the simplest situation taking into account the coupling between 

= and = 1 Landau levels. The strength of the Landau level coupling is parameter- 
ized by 7 = r/hcUc- Our calculated numerical results for the critical exponents(using the 
(5-function potential model) are summarized in Table 4. In Fig. 6 we show one set of represen- 
tative scaling data for the finite 7 Landau level-coupled case. We observe that the calculated 
critical exponents for the lower and the higher energy branches of the same Landau band are 
different due to the strong asymmetry introduced by Landau level couphng effect(£^c 7^ 0). 
This asymmetry is introduced by our model which considers coupling only between the lowest 
two Landau levels. The coupling obviously plays a much stronger role in the higher(lower) 
energy branch of the N — 0(1) Landau band because the couphng is the strongest among 
those states. The most important result of our Landau level coupled calculations is that z/q 
and ui come closer together in the presence of Landau level coupling, allowing us to spec- 
ulate that universality may be restored even for the 5-function scattering potential 
once Landau level coupling is included in the theory. 

We have also carried out calculations for two coupled Landau levels in the presence of 
finite range impurity potentials. These calculations, which are necessarily on somewhat 
smaller system sizes, bear out our basic conclusions that finite potential range and Landau 
level coupling bring uq and z/i close together numerically. We emphasize that by necessity 
our numerical results for the finite range disorder potential and for the Landau level coupling 
case are much less quantitatively accurate than the simple zero-range uncoupled calculations. 
The qualitative trends, however, are reliable and can be trusted. For quantitatively accurate 
results in the presence of finite range disorder potential and Landau level coupling effects, 
one needs to consider substantially larger values of system width which are not currently 
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computationally accessible. We refrain from showing error bars for our finite range disorder 
and Landau level coupling calculations except to mention here that the error bars in Tables 
2-4 are at least an order of magnitude larger than those in Table 1. 

IV. Discussions and Conclusion 

In this paper, we study localization and scaling properties of a disordered two-dimensional 
electron system in the presence of a strong external magnetic field, with emphasis on the 
question of universality of the localization transition in various Landau levels. The impurities 
are treated as random distributed scattering centers with short range 5-function or finite 
range Gaussian potentials. Our localization calculation employs a strip geometry for which 
we obtain the Lyapunov exponent by setting up a suitable transfer matrix. This is done 
rather efficiently in a parallel algorithm. The Lyapunov exponents of the transfer matrix 
are calculated to determine the finite-size localization length. The critical behavior of the 
disordered two-dimensional system near the center of the Landau bands is then studied 
using the standard one-parameter finite-size scaling analysis. We systematically investigate 
the infiuence of varying the impurity concentration, the scattering potential, and the Landau 
level index on the critical exponent. We also study the effect of Landau level couphng on 
the localization exponent. 

In the weak disorder limit, we consider scattering only within a single Landau level, ne- 
glecting any Landau level coupling coupling effect. We find that the critical exponent and 
the scaling property both depend on the Landau level index, particularly for the (5-function 
disorder. Within a single Landau level, especially the lowest Landau level, we observe uni- 
versal scaling behavior independent of impurity concentration, scattering potential type (re- 
pulsive or attractive or both, short range or long range). The calculated critical exponent 
varies a little for different impurity concentrations which we attribute to finite-size or finite- 
concentration effect. We beheve that the dependence of the critical exponent on the Landau 
level index comes from the different symmetry properties of the Landau wavefunction in 
different Landau levels. For short-range scattering potentials with the range smaller than 
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or comparable to the size of the Landau wavefunction, the detailed nature of the Landau 
wavefunction may play a role in the critical behavior, leading to different values of for 
N = and 1. In the presence of only repulsive scattering potential, we observe an expected 
asymmetry in the critical exponent arising from the asymmetry of the density of states which 
is due to the finite impurity concentration effect. For finite range scattering potentials, we 
see the clear trend of decreasing critical exponent with increasing potential range for the 
higher Landau level N — 1 (without much effect on uq). We speculate that for real long 
range scattering potentials with range larger than the size of Landau wavef unctions, the 
dependence of the critical behavior on the details of the symmetry of the Landau wavefunc- 
tion will disappear, producing true universality in the critical localization behavior. This 
percolation limit, beyond the scope of our work, has actually been studied both analytically 
and numerically [21-23]. 

In the strong disorder limit, the Landau level broadening due to the presence of disorder 
is comparable to the inter-Landau-level energy difference, so the coupling among Landau 
levels should be taken into account. It is, in fact, not justified to ignore lower Landau 
levels even in the single- Landau-level approximation when one is considering the localization 
properties of the excited levels(e.5'. A'" = 1). In our numerical calculation, we just consider 
coupling between the lowest two Landau levels. This is, of course, a very crude approximation 
for Landau level coupling which we are forced to make because of our limited numerical 
capability. We find that Landau level coupling, even in this perturbativc sense, has some 
qualitative effect. For the (5-function potential we find that in the strong disorder Landau 
level coupling regime, the numerical values of the critical exponents for different Landau levels 
get closer as we increase the impurity concentration while the scaling properties remain the 
same as in the non-coupling case. The critical energies shift from the center of each Landau 
band because of the asymmetry introduced by Landau level coupling. This new critical 
behavior arises from the mixing of symmetry of different Landau wavef unctions. A similar 
trend of decreasing critical exponent in the higher (A^ = 1) Landau level is seen for finite range 
Gaussian impurities as well. We believe that both finite range disorder and Landau level 
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coupling reduce the role of the detailed symmetry of the Landau wavefunction in the critical 
localization problem and restore universality in the Landau level localization properties. Our 
limited numerical results only indicate a trend and much larger simulations are needed to 
definitely establish this claim. 

There have been several earlier numerical investigations [10-13] of strong-field localiza- 
tion in disordered two-dimensional electron systems. Our parameterization of the random 
disorder potential (as arising from randomly distributed point impurity scattering centers 
characterized by short-range 5-function or finite-range Gaussian potentials) is similar to that 
used by Ando and Aoki [11] who concentrated mainly on the large impurity concentration 
limit (q — 40), neglecting Landau level coupling, and used rather small system sizes (up 
to M = 16 only). Ando and Aoki [11] concluded that i/q — 2 and i/i ~ 4 for 5-function 
scatterers compared with our results — 2.2 ± 0.1 and v\ ~ 5.5 ± 0.5. For finite-range 
scatterers, they found z/q — 2 for d = 2, but did not calculate v\. We also find that z/o(~ 2.2) 
remains unaffected by the finite range of the scattering potential whereas v\ is reduced 
considerably(i/i < 3), bringing it much closer to and thus restoring universality. The 
work of Ando and Aoki was criticized (for having used system sizes far too small to see 
conclusive evidence of scaling) by Huckestein and Kramer [10] who went to larger system 
sizes(M = 64) and concluded that — 2.34±0.04 for short-range white-noise random poten- 
tial. These latter authors [10] did not systematically investigate the exponents in the higher 
Landau levels or the effects of having a finite disorder potential range. They also neglected 
Landau level coupling effects. Our (5-function potential numerical results for the lowest Lan- 
dau level are quantitatively and qualitatively consistent with the findings of Huckestein and 
Kramer. It should be mentioned that Huckestein and Kramer [10] model the disorder po- 
tential somewhat differently from us - instead of working with random distributed impurity 
centers, they use a random white-noise distribution of disorder potential matrix elements 
with a vanishing correlation length. We point out that the largest system sizes(M = 256) 
used in our work are much larger than those used by Huckestein and Kramer. (It is, in 
some sense, essential that we use larger systems because we are interested in investigating 
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higher Landau levels, finite range random potentials, and Landau level coupling effects.) For 
slowly varying random disorder potentials, the strong-field quantum localization problem 
is equivalent to a classical percolation problem, which has been analytically solved [21] to 
produce v — 4/3. If quantum tunneling is included in the percolation calculation [23], one 
gets V — 7/3 ~ 2.33 which is very close to our numerically calculated vq. It is worthwhile 
to mention that in the percolation limit, the Landau levels are necessarily coupled and the 
disorder is necessarily long-range - it is, therefore, reassuring to know that our finite-range 
disorder and Landau level coupled calculations produce uq ~ which are consistent with 
the percolation results. There has also been a numerical calculation [22] near the percolation 
threshold, including effects of quantum interference and tunneling, leading io u = 2.5 ± 0.5. 

Existing numerical and analytical calculations [10-13,21-23] clearly establish that the lo- 
calization critical exponent is around 2.3 independent of the Landau level index for finite- 
range potentials and/or in the presence of Landau level coupling, and for the 5-function 
potential at least in the lowest Landau level. We emphasize that the percolation limit nec- 
essarily includes Landau level coupling and long-range disorder, and is thus complementary 
to the 5-function short-range disorder in the absence of Landau level coupling. The fact that 
both the percolation calculations and the lowest Landau level finite size scaling analysis for 
5-potantials gives the same critical exponent 2.3 provides evidence in favor of universal- 
ity in the Landau level localization. As mentioned before, our finite range and/or Landau 
level coupled numerical calculations are semi-quantitatively consistent with the percolation 
results. On the other hand, there does seem to be a problem with the 5-function random 
potentials in the higher(A'" = 1) Landau levels(in the absence of any Landau level couphng) 
where at least four independent studies [10,11,12,17] (including our own, which employs the 
largest system sizes) find the critical exponent ui to be substantially different from that in 
the lowest Landau level(i/i 5, i^o ~ 2.3). It is, of course, possible (but highly unhkely in 
our opinion) that this lack of universality for the short-range 5-function disorder is purely 
a finite-size phenomenon and future calculations involving larger system sizes will restore 
universality {i.e. ui — uq) even for 5-function disorder. This scenario is unlikely in our view 
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for two reasons: (1) Early calculations involving small system sizes (M = 16) all the way to 
our large system(M = 256, we have some limited results for M = 512) calculations consis- 
tently show ui ^ 2z/o for (5-function disorder (in the absence of Landau level coupling) - the 
quality of scaling ior N — 1 Landau level in our calculations is quite comparable to that for 
N — level, and, therefore, it is unclear how increasing system sizes further would reduce 
i^i by more that a factor of two to bring it in agreement with z/q; (2) for 5-function disorder, 
we see no obvious intrinsic length scale in the problem which demands that we go to system 
sizes substantially larger than this characteristic length scale to obtain the true (rather than 
effective) exponents - thus it is very hard to understand the nature of any crossover behavior 
which may be dominating the finite-size scaling analysis for = 1 producing ui 5 (we do 
expect that one has to use somewhat larger sizes for analyzing the critical behavior of the 
N — 1 level compare to the N — behavior because the free electron Landau wavefunction 
for iV = 1 is more spread out than that for the = level). Based on these considerations 
we feel that (in the absence of Landau level coupling) the range of the disorder poten- 
tial may actually be a relevant perturbation, and the localization transition in higher Landau 
levels may be non-universal for short-range disorder as our numerical calculations indicate. 
We believe that this lack of universality arises from the fact (Fig. 7) that the matrix element of 
the disorder potential projected onto a particular Landau level shows qualitatively different 
spatial behavior for the higher Landau levels compared with the lowest Landau level as a 
function of the range of disorder. This is explicitly shown in Fig. 7. In particular, the spatial 
behavior of the lowest Landau level matrix element remains qualitatively the same (namely, 
it is a monotonically decreasing function of distance) independent of the range of the disorder 
potential (because the wavefunction for the lowest Landau level is nodeless). On the other 
hand, the higher Landau levels being excited states have nodes and, therefore, the spatial 
behavior of disorder matrix elements for the higher Landau levels has a qualitative depen- 
dence on the range of the disorder, showing spatial oscillations for short-range disorder and 
becoming qualitatively similar to the lowest Landau level situation only when the disorder 
range exceeds the spatial extent of the excited wavefunction. Without a real theory for the 
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localization transition, it is difficult to quantify these qualitative considerations. We believe, 
however, that disorder invariably mixes Landau levels and the independent Landau level ap- 
proximation breaks down, making the localization transition universal even for short-range 
disorder as demonstrated by our numerical results in the presence of Landau level coupling 
effects. 

Before concluding we discuss a number of experimental strong-field localization issues 
which are not understood currently. While experiment seems to be in "agreement" with the- 
ory that the localization exponent v ~ 2.3, it is hard to understand why a non-interacting 
theoretical calculation of the localization exponent should agree with the measured experi- 
mental value because, in general, electron-electron interaction is known[24] to be a relevant 
perturbation which changes the universality class. In this respect, it is reassuring that the 
measured localization exponent(«; ~ 0.4) for the fractional quantum Hall transition [25] 
is the same as the integer quantum Hall case. It is important to have some theoretical 
idea about why interaction does not change the universality class in the strong-field situ- 
ation (or, if it does, why the measured value of k is the same for integer and fractional 
situations in apparent agreement with the non-interacting calculation of the localization 
exponent v). The experimentally found dependence[26] of the measured k on the spin de- 
generacy of the Landau level is also not understood at the present time. In particular, k 
for the spin-unpolarized case is found[26] to be half(Ri 0.21) the value of the spin-polarized 
situation(pa 0.43). In a non- interacting localization model, it is unclear how the electron spin 
can be a relevant perturbation. Inclusion of interaction effects, however, leads to a spin de- 
pendence of fermion localization properties in the zero-field case[27]. Whether the observed 
spin dependence of the exponent k is somehow related to interaction effects or not is totally 
unknown at this stage. We suggest a simple scenario for the spin dependence of n which can 
be experimentally tested. Suppose that there is a small (but unresolved) spin-splitting A£' 
in the spin degenerate unpolarized case. In that case, it is possible that, if << ksT, 
the experimental measurement of k measures only an effective k which is smaller than the 
real k because there are two unresolved critical energies separated by a small (but finite) 
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/S.E and one is observing the combined effects of both the locahzation transitions without 
resolving them. (The argument fails if = 0, but the spin-splitting is unlikely to be 
exactly zero.) A simple analysis shows that this scenario leads to a temperature dependent 
effective k which asymptotically approaches the real «;(~ 0.4) in the high temperature limit 
but is lower than the real k in the experimentally feasible regime, approaching, in fact, a 
saturation(z.e.K, = 0) logarithmically in the low temperature limit. Thus, a measurement 
of this effective n in an intermediate temperature range could produce k ~ 0.2 as observed 
experimentally. This proposed scenario should be tested experimentally through detailed 
temperature dependent measurements in the spin degenerate situation. Finally, we mention 
that very recent microwave frequency-dependent measurements [2 8] show clear evidence for 
finite-frequency (/) dynamic scaling of the strong- field localization transition in the integer 
quantum Hall regime. The experimental finding is that for / <; IGHz, a^x peaks broaden 
(at fixed temperature) with increasing frequency, roughly as Ai? ~ p where 7 ~ 0.4(0.2) 
for spin-split (spin-degenerate) peaks. Since one expects 7 = Xjvz where z is the dynamical 
exponent, one concludes that z — \ and v — 2.3(4.6) for spin-split (spin-degenerate) situa- 
tions. The spin-dependence of the microwave experiments is, of course, consistent with that 
of the temperature dependent measurements and remains unexplained (unless our "trivial" 
explanation applies!). The interesting finding of 2; = 1 is consistent with the recent specula- 
tion [2 9] for dirty boson systems in the context of superconductor-insulator transition in thin 
metal films. Statistical transmutation properties of two dimensional quantum systems allow 
mapping between fermion and boson systems, and it is, therefore, tempting to speculate[29] 
that 2; = 1 is a universal feature of all two dimensional quantum phase transitions [30]. One 
problem in this context is that theory [29] predicts a universal value of the critical resistance 
at the metal-insulator transition whereas experimentally [31] the value of the a^x peak seems 
to be non-universal and substantially below the universal theoretical value of jh. One 
should, however, bear in mind that while the theory applies strictly at T = 0, the exper- 
imental measurements are necessarily at finite temperature, and a simple extrapolation to 
T = to obtain the critical conductance may not work. This issue also requires further 
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investigation. 

Our conclusion is that the strong-field two-dimensional Landau level localization is indeed 
universal(ie. vq — z^i, etc.) except for the short-range random disorder potential which, 
in the absence of any Landau level coupling, gives rise to non-universal localization with 
different Landau levels having different localization exponents. Our work indicates, however, 
that inclusion of a finite potential range and/or coupling between Landau levels restores 
localization universality making vq — ui. 
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Table and Figure Caption 



Table 1. Critical exponents for different Landau levels without Landau level coupling 
for weak 5-function impurities. 

Table 2. Critical exponents for weak repulsive 5-function impurities without Landau 
level coupling. 

Table 3. Critical exponents for weak finite range Gaussian impurity-potentials without 
Landau level coupling. 

Table 4. Critical exponents for strong (5-function impurity-potentials considering cou- 
pling between N — and N — 1 Landau levels. 

Figure 1. N = Landau level, for weak 5-function impurities(cj = 8): (a) The renor- 
malized finite system localization length Xm/M as a function of system size M, different 
symbols represent different energies; (b) One-parameter scaling function, with inset: Lo- 
calization length ^{E), v = 2.32 

Figure 2. = 1 Landau level, for weak 5-function impurities(Q = 4): (a) Xm/M vs. M, 
different symbols represent different energies; (b) One-parameter scaling function, with 
inset: Localization length ^(-B), ly — 5.95 

Figure 3. In the presence of weak repulsive 5-function impurities, = 0, Cj = 4, the 
scaling function for the lower(a) and higher(b) energy branches of the Landau band. 
Corresponding localization length in the critical regime are shown in the insets for lower 
energy branch, Ec — — 0.15, — 3.08, and for higher energy branch, Ec — — 0.156, v — 
2.05. 

Figure 4. For weak finite range Gaussian impurity-potentials, N — Landau band, the 
scaling function with localization length in the insets are presented for (a) d/lc — 0.5,Ci — 
8, 1/ = 2.21; (b) d/lc = 1.0, q = 4, i/ = 2.3. 
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Figure 5. For weak finite range Gaussian impurity-potentials, N — 1 Landau band, the 
scaling function with localization length in the insets are presented for (a) d/lc = 0.5, q ~ 
2, i/ = 3.6; (b) d/lc = 1-0, Q = 4, z/ = 2.51. 

Figure 6. For strong 5-function impurities considering coupling between A'^ = and 
N — 1 Landau levels, the scahng function are presented for Ci — 8 (a) N — lower energy 
branch; (h) N — higher energy branch; (c) N — 1 lower energy branch. Corresponding 
localization length in the critical regime are presented in the insets of (a) Ec = 0.029, p — 
1.55; (b) Ec = 0.100, u = 2.24; (c) Ec = 0.872, u = 2.85. 

Figure 7. Shows the spatial behavior of the calculated matrix elements of the impurity 
potential for N — (solid line), N — 1 (dotted hue), and N — 2 (dashed hne) Landau 
levels, and, for four different values of the potential range d/lc = 0(a), 1(b), 2(c), and 4(d). 
The impurity is located at Xi, and lo = \p2^lc where Ic is the magnetic length. 
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Table 2. 
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3.08 
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Table 3. 
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Ec 
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Table 4. 
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